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Erratum; 

The  letter  "v",  used  as  subscript  or  a  superscript, 
should  he  replaced  by  "v"  throughout. 
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SMALL  DEFORMATIONS  SUPERPOSED  ON  LARGE  DEFORMATIONS 
IN  MATERIALS  WITH  FADING  MEMORY* 

by  A.  C.  Pipkin**  and  R.  S.  Rivlin*** 


1*  Introduction 

The  non-linear  stress-deformation  relation  for  materials 
with  memory  has  been  discussed  by  Green  and  Rivlin  [l]  ’  '  ’  ’  and 
Green,  Rivlin  and  Spencer  [2]*  In  the  present  paper  we  derive 
the  restricted  form  applicable  to  small  dynamic  deformations 
superposed  on  large  static  deformations.  Although  this  can  be 
done' by  making  an  appropriate  linearization  of  the  final  result 
of  the  non-linear  theory,  we  avoid  some  of  the  complexities  of 
that  theory  by  introducing  the  linearization  at  an  earlier  stage. 
The  results  obtained  will  apply  in  particular  to  problems  of  wave 
propagation  in  finitely  deformed  viscoelastic  bodies. 

It  is  assumed  that  the  stress  depends  on  the  deformation 
history,  and  that  the  properties  of  the  material  do  not  change 
with  time.  The  stress  components  are  then  hereditary  functionals 
of  the  deformation  gradients.  Because  rigid  rotation  of  a  deform¬ 
ed  body  has  no  effect  on  stress,  other  than  to  rotate  the  stress 
field,  and  since  the  stress-deformation  relation  must  embody 

this  fact,  the  relation  Is  of  a  certain  restricted  form.  This 
form  was  derived  by  Green  and  Rivlin  [l]  by  using  lnvariant- 

**  The  results  presented  in  this  paper  were  obtained  In  the 

course  of  research  sponsored  by  the  Office  of  Naval  Research 
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****  Numbers  in  square  brackets  refer  to  the  References  at  the 
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theoretic  methods,  and  subsequently  Noll  [3]  derived  an  equivalent 
form#  Noll's  derivation  is  simpler  than  that  of  Green  and  Rivlin 
in  that  it  does  not  make  use  of  invariant  theory,  but  the  form 
he  derives  is  not  suitable  for  the  purposes  of  the  present  paper# 
In  Section  3  we  take  the  opportunity  to  present  a  derivation 
which  combines  some  of  the  simplicity  of  Noll's  method  with  the 
explicitness  of  Green  and  Rivlin' s  result# 

When  the  deformation  consists  of  an  infinitesimal  dis¬ 
placement  field  superimposed  on  a  finite  deformation,  the  stress- 
deformation  relation  can  be  linearized  with  respect  to  the 
gradients  of  the  infinitesimal  displacement  field  (Section  4)# 

In  a  material  with  fading  memory,  if  the  finite  deformation  is 
held  fixed  long  enough  so  that  the  original  deforming  process 
has  no  further  effect,  the  non-linear  functionals  specifying 
the  stress  can  be  expressed  in  terms  of  non-linear  functions 
and  linear  functionals  (Sections  5  and  6)  , 

The  functions  appearing  in  the  linearized  formulation 
are  subject  to  restrictions  in  form  if  the  material  has  any 
symmetries  in  its  undeformed  state  (Section  7)*  In  the  case  of 
isotropic  materials (Section  8),  these  restrictions  are  such 
that  the  constitutive  equation  can  be  brought  into  a  simpler 
form  involving  the  classical  strain  tensor  (Section  9)» 
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Consider  a  body  in  its  undeformed  state*  Let  the 
coordinates  of  a  generic  particle  be  X^,  in  a  fixed  Cartesian 
coordinate  system  x.  If  the  body  is  deformed  as  time  progresses, 
then  the  particle  which  was  originally  at  Xi  moves  to  a  new 
position  x^  at  time  t,  and  the  deformation  may  be  specified  by 
giving  x^  as  a  function  of  Xp  and  ts 

xi  =  xi*Xp,1:*  •  (2.1) 

The  derivatives  3x^/8Xp,  called  deformation  gradients,  may  be  used 
as  measures  of  the  deformation  in  the  neighborhood  of  the  particle 

V 

Suppose  that  the  stress  at  a  given  particle  depends  on 
the  deformation  gradients  at  that  particle,  not  only  at  the 
instant  t  considered  but  also  at  all  previous  times*  The  stress 
components  at  time  t  are  then  functionals  of  the  deformation 
gradients  at  the  particle  considered*  To  express  this  idea,  we 
write 

oo 

dij(t)  =  “  T)/3Xqj  t]  *  (2*2) 

*v  =0 

Explicit  dependence  of  the  functionals  on  t  indicates 
that  the  intrinsic  properties  of  the  material  may  change  as  time 
progresses,  whether  the  material  is  subjected  to  a  deformation 
or  not.  This  possibility  will  immediately  be  rejected,  however, 
for  we  asstime  that  the  material  has  the  following  property!  if 
two  deformation  histories  differ  only  by  a  time  shift,  then 
the  corresponding  stress  histories  differ  only  by  that  same  time 
shift*  It  is  easy  to  see,  and  correspondingly  easy  to  show 


3 


Cll-66 

2.  Materials  with  memory  and  hereditary  functional 

Consider  a  body  in  its  undeformed  state*  Let  the 
coordinates  of  a  generic  particle  be  X^,  in  a  fixed  Cartesian 
coordinate  system  x.  If  the  body  is  deformed  as  time  progresses, 
then  the  particle  which  was  originally  at  X.^  moves  to  a  new 
position  x^  at  time  t,  and  the  deformation  may  be  specified  by 
giving  x^  as  a  function  of  Xp  and  ts 

xi  =  xi^Xp>t^  •  (2,1) 

The  derivatives  9x.j/9Xp,  called  deformation  gradients,  may  be  used 
as  measures  of  the  deformation  in  the  neighborhood  of  the  particle 

V 

Suppose  that  the  stress  at  a  given  particle  depends  on 
the  deformation  gradients  at  that  particle,  not  only  at  the 
instant  t  considered  but  also  at  all  previous  times.  The  stress 
components  at  time  t  are  then  functionals  of  the  deformation 
gradients  at  the  particle  considered.  To  express  this  idea,  we 
write 

oo 

d1;j(t)  =  fi;j[dXpCt  t]  •  (2.2) 

Explicit  dependence  of  the  functionals  on  t  indicates 
that  the  intrinsic  properties  of  the  material  may  change  as  time 
progresses,  whether  the  material  is  subjected  to  a  deformation 
or  not.  This  possibility  will  immediately  be  rejected,  however, 
for  we  assume  that  the  material  has  the  following  property*  if 
two  deformation  histories  differ  only  by  a  time  shift,  then 
the  corresponding  stress  histories  differ  only  by  that  same  time 
shift.  It  is  easy  to  see,  and  correspondingly  easy  to  show 


Cll-66 


4 


rigorously*  that  this  assumption  implies  that  the  functionals 
cannot  depend  on  t  explicitly*  Thus  in  place  of  Eq*  (2*2)  we 
write 


su<t)  =  fid[8xp(tT:0x)/aV' 


(2.3) 


The  functionals  are  then  said  to  be  hereditary* 
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3*  SaJaiiaa  sL  J3ia  physical  system 

Not  every  relation  of  the  form  (2*3)  is  physically 
meaningful.  When  an  arbitrary  rigid  rotation  is  superimposed  on 
a  given  deformation  history,  the  stress  field  undergoes  an  equal 
rigid  rotation.  The  form  of  the  constitutive  equation  must  be 
such  that  this  condition  is  satisfied. 

let  be  a  deformation  history  obtained  by  superimposing 
a  rigid  rotation  on  the  deformation  x^t),  so  that 

( JCp,t)  =  aij(t)x;j(Xp,t)  .  (3.1) 

Here  the  coefficients  ^^(t)  must  satisfy  the  orthogonality  condi¬ 
tions 

aik(t)ajk(t)  =  aki(t)ak3(t)  =  bl)  >  '  aij(t)  1#  (3,2) 

From  Eq.  (3*1),  the  deformation  gradients  for  the  new  deformation 
are 

ax1(t)/3Xp  =  a1 j (t )3Xj (t)/axp  .  (3.3) 

Letting  «3 ^ ^ Ct )  denote  the  stress  components  associated  with  the 
deformation  x^Ct) ,  Eq,  (2.3)  implies  that 

°o 

d,  . (t)  =  f^[8x_(t  -  tt)/3X  ]  .  (3#1*-) 

ij  ij  P  ^=0 

The  condition  to  be  satisfied  is  that  the  components 
d"ij(t)  are  related  to  the  components  ^(t)  by  the  rotation  a^^(t)o 
Thus, 


dij(t)  =  ailc(t)a^(t)csk^(t) 


(3.5) 


p 
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By  using  Eqs.  (2.3)  and  (3.4)  in  Eq.  (3.5) ,  and  making  use  of 
Eqs.  (3*2)  and  (3*3) *  we  obtain 

co 


[8VtT:0'0/3X4] 


(3.6) 


oo 


s  aki(tHj(t)fl*t[apr(t  •  'r)axr(t  “  t)/9X^3» 


TssO 


The  functionals  f^  must  be  ofsuahform  that  Eq*(3.6)  is  satisfied 
no  matter  what  the  choice  of  a^Ct)  may  be,  subject  to  the 
restrictions  (3.2). 

We  may  in  particular  let  a^Ct)  be  the  j-component  of 
vector  number  i  in  an  orthogonal  system,  formed  from  the  three 
vectors  QXj/dXj,  3x^/9X2,  and  dx^/dX^  by  the  Schmidt  orthogonali- 
zation  process.  In  the  first  step  of  this  process,  a^(t)  is 
constructed  by  dividing  Sx;J(t)/dX1  by  its  length.  We  thus 
obtain 

alj(t)  =  2  OKjM/oXi  .  (3.7) 


Here  G11(t)  is  one  of  the  set  of  inner  products  Gp(|t),  called 
strain  components,  which  are  defined  by 


W°  - 


9xi(t)  ax^t) 

3X  dX 
p  q 


(3.8) 


Next,  we  construct  the  unit  vector  in  the  plane  of  Qx^/aX^ 

and  3Xj/9X2,  perpendicular  to  3Xj/sX^  and  forming  a  positive 

inner  product  with  dx^/dX^,  Using  the  notation  Ap^(t)  for  the 

cofactor  of  G  (t)  in  the  determinant  G(t)  =  I  G  (t)| ,  this 
pq  P4 

vector  may  be  written 
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a2j(t)  =  [Gn(t)A33(t>]  2[-  G12(t)9x;J(t)/aX1 

+  G11(t)ax;](t)/ax2]  .  (3.9) 

Finally,  a3j(t)  Is  taken  to  be  the  unit  vector  perpendicular  to 
aXj/aX^  and  3x^/aX2  whose  inner  product  with  ax^/aX^  is 
positive.  Thus, 

Sx,(t) 

a^Ct)  =  (t)G(t > ]  *  lA-^Ct)  ■ 

ax,(t)  ax,(t) 

+  A23(t)  ax2  +  ^3(t)  ax3  ^  •  (3,10) 

We  note  that  G^Ct),  d^^Ct),  and  G(t)  are  positive  in  any 
deformation  possible  in  a  real  material.  Equations  (3.7),  (3.9) 
and  (3«10)  may  be  abbreviated  as 

aij(t^  =  Cik^axj^/aXk  i  (3.11) 

where  the  coefficients  Cj^t)  are  factions  of  the  strain 
components  G^Ct). 

Since  the  coefficients  a.jj(t)  constructed  in  the  manner 
described  above  satisfy  the  restrictions  (3.2),  Eq.  (3.6)  must 
be  satisfied  for  this  choice  of  a^Ct).  Thus,  if  Eq.  (3. 11)  is 
used  in  Eq.  (3.6),  and  Eq,  (3,8)  is  used  to  simplify  the  expres¬ 
sion  for  the  functional  argument,  we  obtain 
00 

f,,[axn(t  -  T)/ax J 

13  p  x=0  q 

ax, (t)  ax,(t)  00 

=  ckm(t)  ax m  c<ln(t)  ai  '  fkt^Cpr(t  ~T)Grq(t  “ 

ax, (t)  ax,(t)  00 

°  "5T—  "ST-  pmtopq(t  *  ^  »  <say>* 

“  11  T  =0 


(3.12) 
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The  second  member  of  this  equality  is  of  the  form  shown  in  the 


last  member  because  the  coefficients  C^Ct)  are  functions  of 
the  strain  components  G^(t). 

It  is  thus  seen  that  if  Eq,  (3.6)  is  satisfied,  the 
functionals  f^  can  be  written  in  the  form  (3*12).  It  is  easy 
to  show  that,  conversely,  if  the  functionals  are  of  the  form 
(3.12),  then  Eq.  (3.6)  is  satisfied  for  all  orthogonal  trans¬ 
formations.  The  constitutive  equation  is  therefore  of  the 
restricted  form 

3x* (t)  3x*(t)  °o 

dij(t)  =  ~axT  aFT“  FkjfGpq(t  “  •  (3*13) 

K  *  TT=0 


We  note  that  since  the  stress  tensor  is  symmetric, 
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§aall  deformations  superimposed  on  large  deformations 

Consider  a  deformation  obtained  by  superimposing  a  dis¬ 
placement  field  eu, (X_,t)  on  the  basic  deformation  x.  =  x.(X  ,t): 

1  p  lip 

x*(Xp,t)  =  x1(Xp,t)  +  eu^Xpjt).  (4.1) 


The  stress  components  corresponding  to  this  deformation  are,  from 
Eq.  (3.13), 

ax,  tt;  ax.it;  *  oo 

■ai —  F’^  3 ,  (4.2) 


ax!(t)  axtct) 

su(t>  -  -nj-  “4 


FM.[Gpq(t"*):i’ 

T=0 


where  the  strain  components  G  (t)  are  those  obtained  by  replacing 
* 

x^(t)  by  x^(t)  in  Eq,  (3.8).  We  wish  to  obtain  the  linearized 
form  of  Eq.  (4.2),  valid  for  small  enough  values  of  e  , 

By  using  Eq.  (4.1)  and  neglecting  terms  of  order  e  ,  we 

obtain 


ax*(t)  ax*(t) 
axk  9Xi 


3x1(t)  dX^Ct) 
8Xk  3Xt 


au. (t) 


aXk 
3x^(t) 


axj(t) 


aXk 


ax, 

au.(t) 

— j —  i 

ax,  J 


(4.3) 


ijc 

The  strain  components  Gp^(t)  appropriate  to  the  new  history  can 
be  found  by  setting  i=j  in  Eq.  (4.3).  Thus, 


G^(t)  =  G^t)  +  2eEk<t(t)  . 


(4.4) 


Here  the  functions  G^^t)  are  the  strain  components  for  the  basic 
deformation  xi(t),  and  the  functions  Ek^(t)  are  defined  by 


2Eia(t) 


9u1(t)  ax^(t)  axi (t).  aui(t) 
axk  “I" +  axk  dXl  # 


(4.5) 
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We  now  assume  that  the  functional  is  particularly 

* 

dependent  on  the  value  of  C^Ct-T)  and  its  time  derivatives 

at  t  s  0,  and  is  a  continuous  functional  of  G^Ct-Tr)  over  the 
1  pq 

range  0  <  t  <  °o#  Using  the  notation 

°pq0>  =  Gpq(t)  ’  »  ft*6) 

we  suppose  that  may  be  expressed  as  the  sum  of 

a  number  of  terms,  each  of  which  is  the  product  of  a  differen¬ 
tiable  function  fy,  say,  of  G*^  (y=0,l,2,f •  •),  and  a  continu- 

% 

ous  functional  <J>y,  say,  of  Gpq(t-T?  over  the  range  0  <  x  <  oo  # 
Then, 

PfctC<Spq(t-T>]  *  S  *  (lt*7) 

T=0  V  T=0 

Using  Taylor’s  theorem,  we  have,  with  (4.4)  and  the  neglect  of 
terms  of  order  higher  than  the  first  in  e  , 

f  (G*^)  =  f  (G^)  +  2e  S  df„/8G*^l  ,  (4.8) 

v  pq  vv  pq  J  £  pq  v/uvxpq  |  •  v  * 

In  (4.8)  we  use  the  notation 

Em  =  Epq(t)  >  Epq)  =  *  ft*9) 

By  using  the  integral  representation  for  continuous 

2 

functionals,  it  is  seen  that  neglecting  terms  of  order  8  ,  the 
functional<I>v  may  be  expressed  in  the  form 
*  oo  oo 

°vCGpq(t-T):i  =<tvi:Gpq(t-1?)] 

T=0  oo 

+  e  J<E>^)CGpq^)^3SrS(t-T)drr. 

0  Z=0 


(4.10) 
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By  using  (4,8)  and  (4.10)  in  (4.7),  and  again  neglecting 

o 

terms  of  order  e  ,  we  see  that  Fj^  is  expressible  in  the  form 


00 


CO 


Fl&tCGpq(t-'c):i  =  PfctCGpq(t"‘e)3 


T=0 


T=0 


+e2  aW  [G  Tt-r);  G(fi)]E(v)(t) 
AktrsL  pq  pq  J  rs  v  ' 

T=0 


V 

CO 


+  e  1  \frs[Gp.  J*-55  >  °M>  ’  T3Ers(t-T)aT.  (4.11) 

o  K=o 


The  constitutive  equation  for  small  deformations  super¬ 
imposed  on  large  deformations  is  now  obtained  by  using  Eqs.  (4.3) 

p 

and  (4.11)  in  (4.2),  again  neglecting  terms  of  order  e  •  Since 
e  was  introduced  only  to  facilitate  the  linearization  process, 
it  will  be  omitted  after  this  step;  i.e.  we  take  e =1  with  the 
stipulation  that  the  displacement  gradients  and  their  derivatives 
are  small.  The  linearized  equation  is  then 


(t) 


=  [ 


ax1(t)  ax^t) 

axk  ex^ 


+ 


+ 


au1(t)  ax^t) 

axk  dXi 
axi(t)  au^t) 

axk  ax^ 


CO 


]PUCapq(t-T):l 

T=0 


3x1(t) 

ex,. 


ex,.(t) 

ax  o 


oo 

(t-T) 
pq 
T=0 


G('u)]E(v)(t) 
pq  J  rs  v  ' 


axj (t)  ax,(t) 


»  .'OE -.(t-^dT:  p 
(4.12) 


Since  the  stress  tensor  is  symmetric,  then 

Aktis=AlSs>  and  Kkfcrs=IWs:  Slnee  Ers)=Esr)>  there  is 

loss  of  generality  in  taking  and  Kk/trs=Kk^gr  as  well* 
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If  the  basic  deformation  x^X^t)  were  time -Independent, 
the  functionals  F^j,  and  in  Eq.  (4*12)  would  reduce 
to  ordinary  functions  of  the  constant  strain  components  Gp<^. 

The  basic  deformation  is,  of  course,  necessarily  time-dependent, 


because  it  must  involve  a  transition  from  the  undeformed  state 


to  whatever  fixed  state  is  contemplated*  However,  the  same 
simplification  is  achieved  if  a  material  with  fading  memory  is 
held  in  a  fixed  state  of  deformation  for  so  long  that  the  original 
deforming  process  is  forgotten. 

let  the  basic  deformation  x^(Xp,t)  in  Eq.  (4,1)  be  one 
in  which  the  body  remains  in  a  fixed  state  x^(Xp)  after  time 
t-^j  l.e* 

XfCXpjt)  =  x±(xp),  t  >  o.  (5.1) 


Using  the  notation  G  (t)  for  the  strain  histories,  as  before, 

f  4 

let  denote  the  constant  strains  associated  with  the  fixed 
PQ 

deformation  x1(X  ),  Since  G  (t)  =  G  for  t>0,  then  the 
derivatives  Gp^'  Qx/O)  in  Eq.  (4.12)  all  vanish. 

Vfe  assume  that  after  a  sufficient  amount  of  time  has 


elapsed,  the  deformation  history  before  t=0  has  an  arbitrarily 
small  effect  on  the  stress.  This  assumption  requires  that  the 


following  condition  be  satisfied  by  the  functionals  Fk£  which 

appear  in  Eq.  (4,2): 

oo.  00 

+  -  Fi3[Gp<i+2Epq(t-li]  °* 

T«0  T=0  (5.2) 


As  a  particular  case  of  Eq.  (5.2),  we  have 


As  indicated,  the  functionals  reduce  to  functions  of  the  constant 
values  in  this  case.  By  considering  the  linearized  form  of 
Eq.  (5.2),  i.e.  the  form  corresponding  to  Eq.  (4.12),  it  can  be 
shown  that  the  functionals  and  also  reduce  to 

functions  of  G  at  large  times. 

When  the  deformation  has  been  hold  fixed  long  enough  that 
the  history  before  t=0  has  no  perceptible  effect,  the  functionals 
Fij,  and  ^ijk-t  may  be  rePlaced  by  their  limiting  forms, 

which  are  functions.  The  constitutive  equation  (4.12)  then 
assumes  the  form 

ax,  ax,  au, (t)  ax,  ax,  au,(t) 
dij(t)  =  tsxj  +  axk  ax£  +  ax£  a x^  ^Fk^(GPq) 


ax, 

axk 


ax, 

^i[v 


L(v) 

K'&rs 


(G. 


pq  rs 


OO 

+  I  •  (5*W 

0 

Here  xi  denotes  the  basic  steady  deformation,  with  associated 

strain  components  G„« 

P<1 

In  passing,  we  note  that  if  u,  is  steady,  then  the  combine 
tion  Xj+u^  can  be  regarded  as  a  new  steady  deformation  with  no 
unsteady  part  superimposed.  The  stress  given  by  Eq.  (5*4)  must 
be  the  same,  whichever  way  the  deformation  is  regarded.  From 
this  consideration,  we  find  after  some  manipulation  that  the 
following  relation  must  be  satisfied* 
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°r  9F,  a (G-- ) 

+  I  K13fct(Gpq,,c)dT  “  2 

In  the  remainder  of  this  paper  we  will  consider  only 


those  materials  and  deformations  for  which  Eq,  (5*^)  is  valid. 
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6*  An.-al.ternatlya ..approach  to  materials  with  fading  memory 
In  Sections  2  to  5  we  have  developed  the  constitutive 
equation  for  small  dynamic  deformations  superimposed  on  large 
steady  deformations  in  materials  with  fading  memory#  This  was 
done  by  specializing  the  constitutive  equation  for  materials 
with  an  unspecified  type  of  memory#  An  alternative  approach  is 
possible,  in  which  it  is  assumed  initially  that  the  deformation 
consists  of  a  steady  deformation  x^(Xp)  with  a  displacement  field 
u^X  ,t)  superimposed  on  it#  It  is  then  assumed  that  the  stress 
components  depend  on  the  deformation  gradients  ax^/aX^  and  the 
history  of  the  displacement  gradients  9Up(t)/9Xq,  so  that 

oo 

dij(t)  =  f1;J[aup(t-T)/9Xq  $  9Xp/aXq]  .  (6.1) 

TT=0 

To  find  the  restricted  forms  of  f^  which  will  satisfy 
the  rotation  invariance  principle,  \*e  consider  a  new  total 
deformation  x^Ct)  =  x^  +  u^(t)  related  to  the  original  total 
deformation  x^(t)  =  xi  +  u^(t)  by  a  rigid  rotation,  i#e,  by 
Sq.  (3»1)»  The  stress  components  corresponding  to  the  new 
deformation  are  given  by 

=  fijCaUpCt-Tj/aXq  j  9Xp/9Xq]  .  (6.2) 

The  new  stress  components  5*^ ^ Ct )  must  then  be  related  to  the 
components  d^Ct)  by  the  transformation  (3*5)*  By  using  Eqs# 
(6.1)  and  (6,2)  in  Eq.  (3#5)>  and  making  use  of  Eq#  (3*2),  we 
obtain 

fi;)[eUp(t-T)/9xq  i  axp/axq] 

=  aki(t)a<tJ(t)fk(t[9u  (t-'d/ax  j  a*/ax]  . 

q(6.3) 
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At  this  point,  neither  x^  nor  u^Ct)  have  been  defined* 

Only  their  stun,  x^Ct),  has  been  defined  by  (3*1) •  The  nev  total 
deformation  x^(t)  can  be  decomposed  into  a  steady  part  and  an 
unsteady  part  in  either  of  two  essentially  different  ways: 

(i)  the  steady  part  Xj_  may  be  regarded  as  identical  with  the 
steady  part  x^  of  the  original  deformation,  in  which  case  the 
entire  rigid  motion  is  included  in  the  unsteady  part  u^(t), 

(ii)  the  new  steady  part  x^  may  be  related  to  x^  by  a  time- 
independent  rigid  rotation,  with  u^(t)  related  to  u^(t)  by  the 
same  rotation* 

In  case  (i),  we  take  x^=x^*  The  displacement  field  u^(t) 
is  then  defined  by 

u^t)  =  XjCt)  -  x^  »  a^CtJxjCt)  -  x±  .  (6*lf) 

Equation  (6.3)  then  becomes 

CO 

fjj[aup(t-'u)/dX^  f  aXp/dX  ] 

00 

=  akilt)a^U)flct[aprCl:"*)exr<t",:)'/aXq  '  aVaXq  ’  ax/8V‘ 

T=0 

(6*5) 

By  constructing  the  coefficients  a^Ct)  from  the  gradients 
8xi(t)/aX;},  as  in  Sectioh  3»  we  find  that  f^j  must  be  expressible 
in  the  form 

3Xj  (t)  3x,(t)  *  00 

hi  =  -±r  ~w~  Fk-t[GP,(t-T)  i  vaxq]  • 

*  TT=0 


(6.6) 


C 11-66 


17 

When  Eq.  (6*6)  Is  -used  in  the  basic  invariance  require¬ 
ment  (6,3),  and  Eqs.  (3.1)  and  (3*2)  are  used  to  simplify  the 

$ 

resulting  expression,  it  is  found  that  the  functionals 
must  satisfy  the  condition 

*  00  *  _  co 

FijCGpq(t*^)  *  ax/6Xq]  *  FiJ[Gpq(t“T)  >  dxj/dXql  >  (6*7) 

Tr=0  t=0 

where  the  strain  components  (^(t)  are  defined  by  replacing 
xi(t)  by  x^(t)  in  Eq,  (3.8).  By  making  use  of  Eqs.  (3.1)  and 
(3*2),  it  is  found  that 

5PdU)  =  °P9Ct)  •  (6'8) 

We  now  consider  case  (ii).  The  coefficients  a^  in 
Eqs.  (3.1)  and  (3*2)  are  taken  to  be  constants,  and  the  new 
steady  deformation  is  defined  by 

H  ■  aijx3  •  (6*9) 

By  using  Eqs.  (6.8)  and  (6.9)  in  Eq,  (6.7),  we  find  that  the 

% 

functionals  F^  must  be  of  such  form  that 

*  00  *  00 

Plj[;Gpq(t-',:)  >  8x/8V  “  Pi;jtGpq(t— T)  i  V8X/8V  *  (6>10) 

T=0  T=0 

By  constructing  the  coefficients  a^j  from  the  steady  deformation 
gradients  8Xp/aX^  in  the  manner  described  in  Section  3>  we  find 

l|t 

that  the  functionals  F^  depend  on  the  gradients  ax^/aX^  only 
through  the  strain  components  defined  in  terms  of  them.  By 
using  this  result  in  Eq.  (6.6),  and  the  result  so  obtained  in 
Eq,  (6,1),  we  obtain  finally 
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(6.11) 


9X4 (t)  9x>(t)  OO 

3  IXT  axj""  F^Gpq(t"T)  *  Gpq^  * 

The  displacement  gradients  9ui(t)/8Xj  enter  into  Eq,(6.11) 
through  9x^(t)/aXj  and  Gp^Ct-T),  In  general,  the  displacements 
i^Ct-T;)  cannot  be  small  throughout  the  entire  history  of  deforma¬ 
tion,  since  xi  +  u^(t-iD  =  X^  for  t  sufficiently  large.  Thus, 
before  Eq.  (6,11)  can  be  linearized  with  respect  to  the  displace¬ 
ment  gradients,  it  is 'necessary  to  assume  that  the. material  has  a 
fhdlng  memory  and  that  enough  time  has  elapsed  so  that  the  unde¬ 
formed  state  has  been  forgotten.  With  these  assumptions, 

Eq.  (6.11)  can  be  linearized,  and  the  form  obtained  is  again 
given  by  Eq.  (5*1+)  • 
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Certain  restrictions  must  be  placed  on  the  functions 
and  l*1  if  the  constitutive  equation  is 

to  describe  a  material  which  has  some  symmetry  in  its  undeformed 
state.  let  x  be  a  fixed  rectangular  Cartesian  coordinate  system, 
related  to  the  system  x  by  an  orthogonal  transformation,  so  that 


xi  =  s±fi 

where 


Xi  =  sijXj  >  ui(t)  =  sijuj(t)j 
siks jk  =  skiskj  =  &ij  * 


(7.1) 

(7.2) 


Defining  Gpq,.Epq(t),  and  E^qU)  as  in  Eqs.  (3.8),  (4,5),  and 
(4-. 9)  respectively,  but  in  terms  of  the  new  coordinates  x^  and 
Xi?  and  making  use  of  Eqs.  (7*1)  and  (7.2),  we  obtain 

Gi;j  =  siks^Gk^»  EiJ)(t)  =  siks^E(k^(t)*  (7#3) 

The  components  d^j(t)  of  stress,  measured  with  respect  to  the 
system  x,  are  given  by 

=  siksj.e?!6t(t)  •  (7*4) 

If  a  =  iis^ j II  is  a  symmetry  transformation  for  the 
material,  then  Eq.  (5*4)  remains  valid  if  all  quantities  are 
measured  with  respect  to  the  system  x  instead  of  the  system  x. 

By  using  Eq.  (5*4)  and  the  corresponding  expression  for  d^Ct) 
to  eliminate  <*jj(t)  and  d^^(t)  from  Eq.  (7.4),  a  relation  involv- 
ing  F^j,  ,  and  is  obtained.  From  that  relation  and 

Eqs.  (7*1)  to  (7*3),  it  can  be  shown  after  some  manipulation 
that  the  following  conditions  must  be  satisfied  (for  details 
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of  a  similar  manipulation  in  the  isotropic  case,  see  [1]  or  [2]): 

W  "  >  (7'5) 

Aijrs<5P,)ij&)(t>  -  •tt*Xfi.<0pq)^>(t)  >  <7-6) 

KUrs(Spi>'t®rs(t-'')  =  WjAirsV’^rs^-  (7>7) 

These  relations  must  be  satisfied  identically  for  each  transforma 
tion  s  belonging  to  the  group  of  symmetries  of  the  material.  In 

A/ 

the  following  section  we  will  consider  the  consequences  of 
Eqs.  (7*5)  to  (7.7)  in  the  case  of  isotropic  materials. 
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8.  laajaasls  materials 

The  symmetry  group  for  an  isotropic  material  is  either 
the  full  orthogonal  group  or  the  proper  orthogonal  group, 
depending  on  whether  the  material  has  a  center  of  symmetry  or  not. 
Presence  or  absence  of  a  center  of  symmetry  is  irrelevant  in  the 
present  context,  however,  because  Eqs.  (7*5)  to  (7«7)  are  satis¬ 
fied  identically  if  the  transformation  is  the  central  inversion, 
i.e.  if  a1;j  =  - 

Rivlin  and  Ericksen  [4]  have  studied  the  implications  of 
an  equation  of  the  form  (7.5) •  ^hey  have  shown  that  if  the 
functions  are  single-valued  functions  of  the  variables 
and  if  Eq.  (7*5)  is  satisfied  for  every  orthogonal  transformation 
s,  then  the  matrix  F  =  IlF.  J|  can  be  expressed  in  terms  of  the 

~  **  1 J 

matrix  G  =  IIg^II  in  the  form 


F  =  F.I  +  F,G  +  F  G  , 

~  0«y  ±~  2~ 


(8.1) 


Here  I  is  the  unit  matrix  lib^ll.  The  scalar  coefficients  FgjFj  anl 
F2  are  single-valued  functions  of  three  basic  orthogonal  invari¬ 
ants  of  ft, 

til  S  =  Gii>  tr  J2  =  Gi3Gi:*»  tr  S3  =  GijGjkGki  *  (8,2) 


jje 

If,  in  (7*5),  F^  Is  a  polynomial  in  its  arguments,  then  F  may 
be  expressed  in  the  form  (8.1),  with  FQ,  F^,  and  F2  expressible 
as  polynomials  in  the  invariants  (8a2). 


We  note  that  if  F^  is  a  continuous  function  of  G^,  it  may  be 
expressed  as  a  polynomial  as  closely  as  desired  in  any  closed 
domain. 
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The  limitations  on  Wrs  (t-T)  implied  by  (7*7)  may 
be  obtained  in  more  explicit  form  by  considerations  of  the  type 
used  by  Rivlin  and  Ericksen  [4]  and  by  Rivlin  [£]#  We  use  the 
notation 

Wr*5  =  K1  jr3^Gpq)Ers  •  »*3> 

the  dependence  of  Ers  on  t-T  being  understood.  Then,  Eq.  (7*7) 
may  be  written  in  the  form 

Kij(Spq»®rs)  =  siks^Kl&t(Gpq»Ers)  *  (8#4) 


If  we  assume  that  K^g  is  a  polynomial  in  Gpq,  so  that  is 
a  polynomial  in  Gp(^  and  Eyg,  linear  in  Eps,  it  follows  (Rivlin 
[?])  that  K  =  Hk^ j ||  may  be  expressed  in  the  form 


K  =  K  E  +  K,  (EG  +  GE)  +  ^(EG2  +  G^) 

2 

+  2  K^Oar  EGN)GM  .  (8.5) 

M,N=0  ~~  ~ 

The  coefficients  Kq,  K^,  and  Kj^  are  polynomials  in  the 
invariants  (8.2),  and  functions  of  t  •  If  these  coefficients 
are  regarded  as  single-valued  functions  rather  than  polynomials, 
Eq.  (8 .»+)  is  still  satisfied. 

Since  Eq.  (7*6)  is  of  the  same  form  as  Eq.  (7*7),  its 
implications  are  analogous  to  those  of  (7.7).  We  use  the 
notation 


,E (v)  ) 

*  rs  ' 


(8.6) 


By  an  argument  similar  to  that  used  in  deriving  (8*5),  we  find 
that  A(v)  =  IIaM||  can  be  expressed  in  the  form 
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A(v)  =a(v)e(v)  +A(v)flsCv)04(B(v))  +A(v)(e(v)g2^,2e(v)) 

O  /V  X  <v  ej  ~~  2  M  «w  ou  /V 


+  s  A^Ctr  e(v)GN)GM  . 
M.N=0  m  ~  ~  ~ 


(8.7) 


The  coefficients  A^v) ,  a[v)  ,  A^v) ,  and  A^  are  polynomials  in 
the  invariants  (8,2).  Again*  the  invariance  requirements  are 
still  satisfied  if  these  coefficients  are  regarded  as  single¬ 
valued  functions  rather  than  polynomials. 

When  Eq.  (8.5)  is  used  in  (8.3),  (8.7)  is  used  in  (8.6), 
and  (8.1),  (8.3),  and  (8.6)  are  then  used  in  Eq.  (5*4-),  the 
constitutive  equation  for  isotropic  materials  is  obtained.  It 
may  be  VTitten  as 

dij(t) 

ax.  ax.  au.(t)  ax.  ax.  au.(t) 

=  [a\  aT^  +  -8*7-  <  +  ^r]CFo6kt +  FiGk*  +  F2GkpW 

5Xj  /  \  /  \ 

+  axj  aij  2  [a0  +  a^  (6^  +  6tqGkp) 


OO 


+  A2V,^kp°-txGrq  +  %°krGrp^ 

H,N=0  ^  ~  ~  PQ  pq 


ax.  ax 


+  *-*  1  tK« 


axk  ax^  LVT)6kp6<tq  +  Kl(T)(5kp%  +6^qGkp) 


2 

+  z 

H 


^  <8-8> 


Dependence  of  the  scalar  coefficients  on  tr  GN  (N=1.2.3)  is 

«v  '  * 

understood. 
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9.  Afl-jatazmfrlY§  Xgra  to  agnajaiaSlga  afluatton.  .to  iaatoBls 

aa.to.rial? 

Equation  (8.8)  can  be  written  in  a  slightly  simpler  form. 

As  the  first  step,  the  displacement  gradients  are  written  in 
terms  of  spatial  derivatives,  thus 


au^  C  t ) 
axi 


au^t)  ax. 

* 


(9.1) 


Here  the  derivatives  9x^/8Xj  are  deformation  gradients  for  the 
basic  steady  deformation.  By  using  (9.1)  in  the  definition  (4.5) 
of  E  ,  we  obtain 


(9.2) 


where  e^  is  the  classical  strain  tensor  defined  by 


•»  3Uj  (t)  9u,(t) 

eijCt)  “  2^  dxi  +  dx±  1  * 

(9.3) 

The  time -derivatives  E^^  are,  from  (9.2), 

(9.4) 

where 

ei (t)  =  dH'ei;j(t)/dt,A  . 

(9.5) 

In  the  definition  (4.9)  of  ,  and  now  in  Eq, 

differentiation  d/dt  holds  fixed.  Thus 

(9.5),  the 

.1  =  -SL  +  JL  . 

dt  at  at  ax1  • 

(9.6) 

When  this  expression  is  used  in  (9.5) ,  terms  arise  which  are 
quadratic  in  the  displacements.  These  will  be  neglected.  Thus, 
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in  (9*5) >  d/dt  may  be  regarded  as  the  partial  derivative  with 
x^  held  fixed* 

Next,  consider  the  tensor  g  =  llg^ll,  defined  by 


9x4  Qx 


gi3  =  »Xp  «p  * 


(9.7) 


By  using  this  definition  with  Eqs.  (9.2),  (9.^) »  and  the  defini¬ 
tion  (3.8)  of  ,  we  obtain  the  relations 
pq’ 

to  0H  =  to  gH  ,  to(E(v)GN)  =  to(e(v)giW) ,  (9.8) 

M  rst  *  M  «v 


^  (~N)^  =  (|N+1) 


±i 


^  ax^E(v)£N)fct  =  (I£(v)|N+1)i3  » 


(9.9) 


(9.10) 


ax.  ax.  au4 (t)  ax4  ax.  au.(t)  N 

ca5^  sq+  a5^  +  35£  Jk-t 

■  t'W’j.T  ^,k(tW*t+ 

In  the  last  equation,  the  notation  au^/dx^  =  bas  been  used. 

By  making  use  of  Eqs.  (9.8)  to  (9.11) »  the  constitutive 
equation  (8.8)  can  be  written  in  the  form 

sij(t)  =  bik'-ji-  ui,k(t)tH+  uj,i(t)6i^(Fol  +  Fi£2  +  Fs£3)ki 
+  2[A<v)g9(T)g  ♦  A<v)(ge(v)g2  +  g2e(v)g) 

••  O  MN  M  IVfV  A*  ^  ^  ++ 


2 
S 

M.N=0 


A<v)(ge(v)g3  +  g3e(v)g)  _  A^>tr(e(T)gN+1)gM+1] 


'ij 


OO  “I"  w 

f[K  (T)geg  +  K,(T)[geg2+g2eg)  +  K?(t)  (geg3+g3eg) 


(9.12) 


11,N=0 
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The  time-dependence  of  e  and  e^is  understood.  According  to 

n  »v 

(9*8),  the  scalar  coefficients  in  (9*12)  may  "be  regarded  as 
functions  of  tr  gN  (N=l,2,3)  rather  than  of  tr  G^1  (N=l,2,3)» 

The  integration  with  respect  to  time  in  (9.12)  is 
integration  with  held  fixed.  However,  with  neglect  of  terms 
quadratic  in  the  displacements,  it  may  be  regarded  as  an  inte¬ 
gration  holding  x^  fixed. 

•5 

By  using  the  Cayley-Hamilton  theorem,  gJ  may  be  expressed 

A, t 

2 

as  a  linear  combination  of  g  ,  g,  and  I,  with  coefficients  which 

A*  r>J  ^ 

N  / 

are  polynomials  in  the  invariants  tr  g  (N=l,2,3).  Therefore 

r>j 

we  may  write 


Fo£  +  Pl|2  +  F2g3  =  +  fll  +  f2|2  > 


(9.13) 


where  fo,  f^,  and  f2  are  functions  of  tr  gN  (11=1,2,3). 

Denoting  the  integrand  in  Eq,  (9.12)  by  and  its  matrix 
by  k  =  l!k.  .  !l ,  we  see  that  k  is  a  symmetric  matrix  polynomial  in 

A*  lj  A* 

the  symmetric  matrices  g  and  e,  linear  in  e.  The  scalar 

/s*  Aw  a/ 

coefficients  are,  of  course,  not  necessarily  polynomials. 
According  to  Rivlin  [?],  k  may  therefore  be  expressed  in  a  form 
similar  to  (8.5).  In  that  expression  we  replace  K  by  k,  G  by  g 

A#  (V  N  A* 

and  E  by  e  to  get  the  desired  result.  Similarly,  the  summand  in 


fV  Arf 

(9.12)  can  be  expressed  in  a  form  analogous  to  Eq.  (8.7),  by 
replacing  G  by  g  and  E  by  e  in  the  latter  expression.  These 

A*  a*  Arf 

results,  taken  with  Eq.(9,l3),  imply  that  Eq,  (9.12)  can  be 
reduced  to  the  following  form: 
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2- 


=  [&ik£>.j,t+  ui,k^6j-t  +  uj,^t^6ik^fo£  +  fl£  +  f2&  ^k£ 

+  St^^ik6;^  +  alV)(6iksH+6Hgik) 

+  a2  ^ikg  jngri't+&j'tsingnk^ 


M,l=0  aMN)(gM)ij(£N)kt]ek\)(t) 


oo 


+  J  ^0(T)6lk6j/t  +  ^(t)  (6ikgj't+6j£gik) 

0 

+  k2U)(6ikgjngnt+5jtgingnk) 

+  m|=o  U(*\l]ekl(t-T)dr- 

The  scalar  coefficients  are  functions  of  tr  gN  (N=l,2,3)« 


(9.1^) 
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